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Abstract 

We present a class of solutions to the discrete Painleve-II equation for particular 
values of its parameters. It is shown that these solutions can be expressed in terms of 
Casorati determinants whose entries are discrete Airy functions. The analogy between 
the r function for the discrete Pn and the that of the discrete Toda molecule equation is 
pointed out. 
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1. Introduction 



The six Painleve transcendents are of very common occurence in the theory of in- 
tegrable systemsfl]. Nonlinear evolution equations, integrable through inverse scattering 
techniques, were shown to possess one-dimensional (similarity) reductions that are just 
Painleve equations. This feature of integrable PDE's eventually evolved into an integrablity 
criterion[2], the Painleve property being intimately linked to integrablity. Discrete inte- 
grable systems have recently become the focus of interest and an active domain of research. 
The study of partition function in a 2-D model quantum gravity [3,4] led to the discovery 
of the discrete analogue of the Painleve- 1 (Pi) equation. It was followed closely afterwards 
by the derivation of the discrete Pn in both a quantum gravity setting [5] and as a simi- 
larity reduction of a lattice version of mKdV equation[6]. The remaining discrete Painleve 
equations (dPni to dPy) were derived[7] using a more direct approach reminiscent of the 
Painleve- Gambier [8] method for the continuous ones. This method, derived in [9] and 
dubbed singularity confinement, is the discrete equivalent of the Painleve approach and 
offers an algorithmistic criterion for discrete integrability. One important result of these 
investigation is that the form of the discrete Painleve equations are not unique: there exists 
several discrete analogues for each continuous Painleve equations. 

The continuous Painleve equations were shown to be transcendental in the sense that 
their general solution cannot be expressed in terms of elementary functions [10]. In fact, 
this solution can be obtained only through inverse scattering methods. However in some 
particular cases (for special values of parameters) the solution to the Painleve equations 
can be expressed in terms of special functions[ll,12,13]. For example, Pn 

w xx - 2w 3 + 2xw + a = , (1) 

has a solution for a = — (2iV + 1) 

»=0og^),, (2) 
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where tn is given by an iV x N Wronskian of the Airy function, 



Tat = 



Ai 



d 



Ai 



dx 

dx y dx' 



Note that Ai is the Airy function, satisfying 

d 2 



d 



(—) 

^dx' 



Ai 



(if * 



d 



Ai 



(3) 



dx' 



Ai = xAi 



(4) 



From the close analogy that is known to exist between the continuous and discrete Painleve 
equations one would expect special function-like solutions to exist for the discrete Painleve 
equations as well. This is indeed the case. As was shown in [14], dPn has elementary 
solutions that can be expressed in terms of the discrete equivalent to the Airy function. 
In [14] only the simplest of these solutions was derived explicitly. The method for the 
obtention of the higher ones was based on the existence of an auto-Backlund transform for 
dPn but it is not clear how one can get the general expression of these special function 
solutions following this method. In this letter we intend to present the answer to this 
problem. Using Hirota's bilinear formalism we can show that these particular solutions to 
dPn can be written as Casorati determinants whose entries are the discrete analogues of 
the Airy function. 

2. Special Solutions of dP n 



We consider dPn, 



w n+1 + ty n _i = 



(an + P)w n + 7 
l-wl 



(5) 



where a, (3 and 7 are arbitrary constants. First, let us seek a simple solution of eq.(5). It 
is easily shown that if w n satisfies the following Riccati-type equation, 

w n - (an + b) 



w n+ i 



l+w n 



(6) 



then it gives a solution of eq.(5) with the constraint 7 = —a/2. In fact, we have from 
eq.(6) 

w n + (an -a + b) 
w n -i = ; . (7) 



Adding eqs.(6) and (7), we obtain 



w n+1 + W n -i 



l-w n 



(2an — a + 2b + 2)w n — a 



1-wt 



(8) 



which is a special case of eq.(5). Now we put 



w. 



G r 



(9) 



Substituting eq.(9) into eq.(6) and assuming that the numerators and the denominators of 
both sides of eq.(6) to be equal, respectively, we have 



F n+1 = F n - (an + b)G n , 



(10a) 



Gn+l — G n + F n . 
Eliminating F n from eqs.(lOa) and (10b), we see that G n satisfies 



(10b) 



G n +2 - 2G n+ i +G n = -(an + b)G n , 



(11) 



which is considered to be the discrete version of eq.(4) and has a solution given by the 
discrete analogue of the Airy function. By means of the solution, w n is expressed as 



G 



Wr 



n+1 



G n 



- 1 



(12) 



It is possible to construct a series of solutions expressed by the discrete analogue of 
the Airy function. We here give the result, leaving the derivation in the next section. We 
consider the r function, 



A 



n+l 



A n +3 



A n +N-1 A-n+N+1 
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A n +2N-1 



ln+3AT-3 



(13) 



where A n satisfies 



A n+2 = 2A n+1 - (pn + q)A n 



(14) 



We can show that satisfies the following bilinear forms, 



and 



n-l n+2 _ n-l n+2 _ n n+l /-. r\ 

N+l N—l ~ N N 'N'n > \ 10 ) 



r n N \\r n N +1 ~ 2 r^r^ 2 + (pn + ? ) r^ +1 r-+ 3 = , (16) 



^+1^ = -(p(n + 2N) + q) r^+V^ 1 + (pn + q) . (17) 



Applying the dependent variable transformation as 



n+l n 

N+l N , /10 n 

Wn = +1 " 1 , (18) 

r N+l r N 



we obtain a special case of dPn, 



(2pn + (2N - l)p + 2g) w n - (2N + l)p 

w n+ i +w n _i = 2 • ( 19 ) 

1 - to. 



'n 



We note that eq.(8) and its solution is recovered by putting p = a, q = b + 1, and iV = 0. 
We also note that eq.(19) reduces to eq.(l) with a = —(2N + 1) if we choose p = — e 3 , 
q = 1, w n = ew and n = =|, and take the limit of e — > 1. 



3. Derivation of the results 



In this section we show that eq.(13) really gives the solution of eq.(19) through the 
dependent variable transformation (18). 

First, let us prove that the r function (13) satisfies the bilinear forms (15)-(17). For the 
purpose we show that eqs.(15)-(17) reduce to the Jacobi identity or the Pliicker relations. 
Before doing so, we give a brief explanation on the the Jacobi identity. Let D be some 
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determinant, and D ^ J be the determinant with the i-th row and the j-th column 
removed from D. Then the Jacobi identity is given by 



Dl ') dI*) !)(]) D (k 



D D 



i k 



(20) 



J J V ' J W \J J \3 1 

It is easily seen that eq.(15) is nothing but the Jacobi identity. In fact, taking t^_\ as D, 
and putting i = j = 1, k = I = N + 1, we find that eq.(15) reduces to eq.(20). Hence it is 
shown that eq.(13) satisfies eq.(15). 

Let us next prove eq.(16). Notice that rjy is rewritten as 



T N — 



A n 


■ ■ A n+2 N 




■ ■ A n+2 N 


An+N-1 ■ 


■ ■ A n+3 N 


A n 


■ ■ A n+2 N 


An+l 


■ ■ A n+2 N 


An+N-1 ■ 


■ ■ A n+3 N 



A n 
A n +i 



2A 



n+2N-3 



2A n+ i 
2A n+2 - pA n+1 



A n+N _i 2A n+N — (N — l)pA n+N _i 

A n +1 ■ ■ ■ A n+ 2N-3 



2A n+ 2N-3 
2A n+2 N-2 — pA n+ 2N-3 

2A n+3A r_ 4 — (N — l)pA n+3N -5 



= 2 



JV-1 



B (o) 



R^" 1 ) A A 



where B n k \ k = 0, 1, • • •, are given by 



(21) 



Si 0) = A, 



= A n+k + k 4B^ for k > 1 



2 n 



(22) 



Similarly, we have 



(pn + q) T % = 2 



N-l 



A n +i 

An+2 



B 



(o) 

n+2 

o(l) 
n n+2 



(N-l) 



An+3 
Ari+4 

{■n+N+2 



A n +2N-3 
A n +2N-2 



■■ A 



n+3N-4 



(23) 
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Let us introduce the notations as 



3 = 



( A n+j \ 
An+j+l 



V" = 



/s (0 2\ 

n+j 



= 



/ \ 



Vi/ 



(24) 



V / 

For example, and (pn + ?)r]^ are rewritten by 

t]!J- = |0, 2, • • • , 2N — 2| = |0, 2, • • • , 2N — 2, 0| 
= 2 Ar " 1 |0 / , 1,3, ••-,2AT-3| , 
(pn + Q)TN = 2 7V - 1 |1,2 , ,3,5,---,2A^-3| = 2 N ~ 1 \1, 2', 3, • • • , 2N - 3, <f>\ , 
respectively. Now consider the following identity of (2N + 2) x (2N + 2) determinant, 
-1 0' i 1 ••• 2N -5 i i 2iV-3 (p 

-1 0' i 1 ••• 2N -5 i 2iV-3 6 



= 



(25) 



Applying the Laplace expansion on the right hand side of eq.(25), we obtain 

=| - 1, 0', 1, • • • , 2N - 5| x |1, • • • , 2N - 5, 2N - 3, <f>\ 
-| — 1, 1, • • • , 2JV — 5, 2JV — 3| x |O',l,---2iV-5,0| 
+ | - 1, 1, • • • , 2N - 5, <j>\ x |0', 2N — 5, 2N — 3| , 



(26) 



which is nothing but the special case of the Pliicker relations. Equation (26) is rewritten 
by using eqs.(21) and (23) as 



= (p(n -2) + q) r"- 2 t% + _\ - 2 r^ 1 r n N _, + r n N r n N ~_\ , 



(27) 



which is essentially the same as eq.(16). 

We next prove that eq.(17) holds. We have the following equation similar to eqs.(21) 
and (23); 



(p(n + 2N) + q)r™ +2 = -|2, • • • , 2N - 2, 2N + 2\ + 2 7V " 1 |2', 3, • • • , 2N - 3, 2N + 1| . (28) 



Then the right hand side of eq.(17) is rewritten as 



= 



|2, •••,2iV-2,2iV + 2| x |1,3, • • • , 2iV - 1| 
-2 7V " 1 |2',3, •••,2AT-3,2AT+1| x |1, 3, • • • , 2N — 1| 

+2 Ar - 1 |l,2 / ,3,---,2iV-5,2iV-3| x |3, 5, • • -2N — 1, 2N + 1| . (29) 

From the identity 

1 2' i 3 • • • 2N — 3 i i 2N - 1 2iV + 1 

1 2' | j 3 • • • 2iV - 3 i 2N - 1 2N + 1 

= |1,2', 3, - • -,2JV — 3| x |3, 5, • • -2iV — 5,2JV — l,2iV + l| 

- |1, 3, • • • 2JV — 3, 2JV — 1| x |2',3, • • • , 2iV - 3, 2N + 1| 

+ |1,3,---,2AT-3,2AT + 1| x |2', 3, • • • 2AT - 3, 2N - 1| , (30) 

the second and third terms of eq.(29) yield 

— 2 JV_1 11, 3, •••,2AT-3,2AT + 1| x |2', 3, • • • , 2N — 3, 2N — 1| 
= - |1, 3, • • • , 2N - 3, 2AT + 1| x |2, 4, • • • , 2N - 2, 2N\ . (31) 

Hence, eq.(17) is reduced to 

|2,4, •••,2AT-2,2AT + 2| x |1, 3, • • • , 2N — 3, 2N — 1| 
-|1,3,---,2AT-3,2AT+1| x |2, 4, • • • , 2N - 2, 2iV| 

=|1, 3, - - - , 2iV — 1, 2AA + 1| x |2, 4, • - • , 2JV — 2| , (32) 

which is again nothing but the Jacobi identity (20). In fact, taking D = |1, 3, • • • , 2N — 
1, 2N + 1|, i = 1, j = N + 1, k = N and I = N + 1, we see that eq.(20) is the same 
as eq.(32). This completes the proof that the r function (13) satisfies the bilinear forms 
(15)-(17). 
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Finally, let us derive eq.(19) from the bilinear forms (15)-(17). We introduce the 
dependent variables by 



n _n n+3 

n _ T N+1 n _ T N T N /ooN 

V N - n > N n+1 n+2 ' l dC V 
r 7V r 7V r 7V 



Then eqs. (15)-(17) are rewritten as 



Vn' 1 = <-\ ( 1 - 4=1 ) > (34) 



n+2 o„,n+l 
N 



2< 1 + (pn + g )u>; = 0, (35) 



= vp~_\ (-(p(n + 2 AT) + g) + (pn + ^ , 
respectively. Eliminating w^r and vn-i from eqs.(34)-(36) and introducing w n by 



(36) 



„n+l 

= -Sr- - 1 , (37) 



we obtain eq.(19). 



4. Concluding Remarks 



In this letter, we have discussed the solution of dPn, (for semi-integer values of the 
parameter 7/ct ) and shown that it is expressed as a Casorati determinant of the discrete 
Airy function. The most remarkable result is the structure of the r function (13). The 
subscript of A n does not vary in the same way in the horizontal and vertical directions: it 
increases by one with each new row and by two with each new column. This is a feature 
which has not been encountered before in other discrete integrable systems. 

Before concluding let us point out the relation with the Toda molecule equation. It 
is known in genaral that the r function of Pn satisfies the Toda molecule equationfll], 

d 2 d 2 

—2 T N ■ T N - (— T N ) = Tat+iTtv-i , N = 0, 1, 2, • • • , (38) 
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whose solution is expressed as 



tn 



f d f 

JL f i-^-S 2 f 



( d \ N ~ X f 

( d ) N f 



( d O^- 1 / id \n f 



( d \2iV-2 



(39) 



where / is an arbitrary function. It is clear that eq.(3) is a special case of eq.(39). Hence, 
we may expect that the r function of dPn satisfies the discrete Toda molecule equation 
proposed by Hirota[15], 



A 2 t n ■ t n - (At n ) 2 = t n+1 t^ + _\ , N = 0,1,2 ■ 



(40a) 



or 



n+2 n _ I n+l\ 2 _ n n+2 at _ n i n 

T N T N \'N ) — T N+1 T N-1 > JV — V> L > Z , 



whose solution is given by 



TN = 



fn A/ n 

A /„ A 2 /„ 



A^" 1 /„ 



N-l 



fn A N f, 



.2N-2 



(40b) 



(41) 



for arbitrary f n , where A is a forward difference operator in n defined by 

n. r N — r N — r N . 

However, because of the difference of the structure of the r function mentioned above, that 
of dPn does not satisfy the discrete Toda molecule equation (40) itself. In fact, eq.(15) 
may be regarded as an alternative of eq.(40), which is rewritten as 



(AAV#) • r£ - (Ar N ) ■ (A'r N ) = r N+1 r^_\ 



(43) 



where A' is given by 



A' -r n — n- n + 2 n~ n 
N ' JV ~~ T N 
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Indeed, eq.(43) also reduces to the ordinary Toda molecule equation (38) in the continuum 
limit. 

It is expected that the other discrete Painleve equations have also the solutions ex- 
pressed by Casorati determinants whose entries are the discrete special functions. In 
particular for dPni it was shown in [16] that solutions in terms of discrete Bessel functions 
exist for some values of the parameters, while for dPiv the particular solutions are in terms 
of discrete parabolic cylinder (Weber- Hermite) functions [17]. In a forthcoming paper we 
intend to present Casorati determinant-type solutions for these discrete Painleve equations. 
One more interesting point concerns the existence of rational solutions. Both continuous 
and discrete Painleve equations possess such solutions and it should be in principle to 
obtain general expressions for them in terms of Casorati determinants. 

One of the author(K K) was supported by Grant-in- Aids for Scientific Research Fund 
from the Ministry of Education, Science and Culture of Japan (3515). 

References 

[1] Painleve Transcendents 1992, Ed Levi D and P Winternitz, NATO ASI series B 

(New York: Plenum) 
[2] Ablowitz M J, Ramani A and Segur H 1978, Lett. Nuov. Cim. 23 333 
[3] Brezin E and Kazakov V A 1990, Phys. Lett. 236B 144 
[4] Its A R, Kitaev A V and Fokas A S 1990, Usp. Mat. Nauk 45,6 135 
[5] Periwal V and Shevitz D 1990, Phys. Rev. Lett. 64 135 
[6] Nijhoff F and Papageorgiou V G (1991), Phys. Lett. A 153 337 
[7] Ramani A, Grammaticos B and Hietarinta J 1991, Phys. Rev. Lett. 67 1829 
[8] Ince E L 1956, Ordinary Differencical Equations (New York: Dover) 
[9] Grammaticos B, Ramani A and Papageorgiou V G 1991, Phys. Rev. Lett. 67 1825 
[10] Painleve P 1902, Acta. Math. 25 1 

12 



[11] Okamoto K 1987, Ann. Mat. Pura Appl. 146 337; 1987, Japan J. Math. 13 47; 1986, 

Math. Ann. 275 221; 1987, Funkcial. Ekvac. 30 305 
[12] Nakamura A 1992, J. Phys. Soc. Jpn. 61 3007 

[13] Gromak V A and Lukashevich N A 1990, The Analytic Solutions of the Painleve 

Equations (Minsk: Universiteskoye Publishers), in Russian 
[14] Ramani A and Grammaticos B 1992, J. Phys. 25A L633 
[15] Hirota R 1987, J. Phys. Soc. Jpn. 56 4285 

[16] Grammaticos B, Nijhoff F, Papageorgiou V G, Ramani A and Satsuma J, submitted 
to Phys. Lett. A 

[17] Tamizhmanni K M, Ramani A and Grammaticos B, to appear in Lett. Math. Phys. 



13 



